Search results for "Gaussian quadrature"
showing 5 items of 5 documents
Quadrature effects on the accuracy of flux calculations in realistic atmospheres
1993
Abstract We have investigated the accuracy of five different quadrature methods—equal steps in θ, equal steps in cos θ, Gaussian, double Gaussian and Gauss-Lobatto—on the accuracy of fluxes in realistic aerosol atmospheres, using the Gauss-Seidel method. In addition, a range of Gaussian quadrature stream numbers from two to 32 were compared. The atmospheric models considered are those recently presented by Lenoble, with the exception that we have used Henyey-Greenstein phase functions in place of Mie. Our results should be easily reproduceable by any other workers interested in similar realistic atmospheres. A table of Gauss-Lobatto weights and points is provided as an appendix.
Error Bounds for the Numerical Evaluation of Integrals with Weights
1988
This paper is concerned with a procedure of obtaining error bounds for numerically evaluated integrals with weights. If \( - \infty \mathop < \limits_ = a < b\mathop < \limits_ = \infty \), w integrable over [a,b] and positive almost everywhere, then an approximation of \({I_W}f: = \int\limits_a^b {w\left( t \right)f\left( t \right)dt} \) by a quadrature rule \({Q_n}f: = \sum\limits_{i = 0}^n {{\alpha _i}f\left( {{t_i}} \right)} \) is leading to the error Enf ≔ Iwf ‒ Qnf. An algorithm is derived for the computation of bounds for |Enf| depending on the smoothness of the integrand f and on the degree of exactness of Q. As initial values this algorithm needs moments of the weighting function w…
Quadrature Formula Based on Interpolating Polynomials: Algorithmic and Computational Aspects
2007
The aim of this article is to obtain a quadrature formula for functions in several variables and to analyze the algorithmic and computational aspects of this formula. The known information about the integrand is {λi(f)}i=1n, where λi are linearly independent linear functionals. We find a form of the coefficients of the quadrature formula which can be easy used in numerical calculations. The main algorithm we use in order to obtain the coefficients and the remainder of the quadrature formula is based on the Gauss elimination by segments method. We obtain an expression for the exactness degree of the quadrature formula. Finally, we analyze some computational aspects of the algorithm in the pa…
Time-harmonic elasticity with controllability and higher-order discretization methods
2008
The time-harmonic solution of the linear elastic wave equation is needed for a variety of applications. The typical procedure for solving the time-harmonic elastic wave equation leads to difficulties solving large-scale indefinite linear systems. To avoid these difficulties, we consider the original time dependent equation with a method based on an exact controllability formulation. The main idea of this approach is to find initial conditions such that after one time-period, the solution and its time derivative coincide with the initial conditions.The wave equation is discretized in the space domain with spectral elements. The degrees of freedom associated with the basis functions are situa…
Gaussian quadrature rule for arbitrary weight function and interval
2005
Abstract A program for calculating abscissas and weights of Gaussian quadrature rules for arbitrary weight functions and intervals is reported. The program is written in Mathematica. The only requirement is that the moments of the weight function can be evaluated analytically in Mathematica. The result is a FORTRAN subroutine ready to be utilized for quadrature. Program summary Title of program: AWGQ Catalogue identifier:ADVB Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADVB Program obtained from: CPC Program Library, Queens University, Belfast, N. Ireland Computer for which the program is designed and others on which it has been tested: Computers: Pentium IV 1.7 GHz processor Ins…